Influences of degree inhomogeneity on average path length and random walks in 

disassortative scale-free networks 
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Various real-life networks exhibit degree correlations and heterogeneous structure, with the latter 
being characterized by power-law degree distribution P(k) ~ fc -7 , where the degree exponent 7 
describes the extent of heterogeneity. In this paper, we study analytically the average path length 
(APL) of and random walks (RWs) on a family of deterministic networks, recursive scale-free trees 
(RSFTs), with negative degree correlations and various 7 G (2, 1 + jjrf], with an aim to explore 
the impacts of structure heterogeneity on the APL and RWs. We show that the degree exponent 
7 has no effect on the APL d of RSFTs: In the full range of 7, d behaves as a logarithmic scaling 
with the number of network nodes N (i.e. d ~ In TV), which is in sharp contrast to the well-known 
double logarithmic scaling (d ~ In In TV) previously obtained for uncorrelated scale-free networks 
with 2 < 7 < 3. In addition, we present that some scaling efficiency exponents of random walks are 
reliant on the degree exponent 7. 

PACS numbers: 89.75.Hc, 89.75.Fb, 05.40.Fb 



I. INTRODUCTION 



The last decade has witnessed tremendous activities 
devoted to the characterization and understanding of 
real-life systems in nature and society [l], H, H, 0| ■ Ex- 
tensive empirical studies have revealed that most real 
networked systems exhibit scale- free behavior Q, which 
means that these systems follow a power-law degree dis- 
tribution P(k) ~ fc -7 with degree exponent 7 6 [2,3]. 
Generally, we call a network with scale-free behavior 
a scale-free network (SFN), which has a heterogeneous 
structure encoded in the characteristic degree exponent 
7: the smaller the 7, the stronger the heterogeneity of 
the network structure. The inhomogeneous degree dis- 
tribution of a SFN has a profound effect on almost all 
other aspect of the network structure. For example, it 
has been established that scale-free behavior is relevant 
to average path length (APL) Q in uncorrelated random 
SFNs, i.e., the APL d(N) for a network with node num- 
ber N depends on 7 when 7 = 3, d(N) ~ ln/V; 
when 2 < 7 < 3, d(N) ~ lnlniV. 

As known to us all, the ultimate goal of studying net- 
work structure (e.g. degree distribution) is to understand 
how the dynamical behaviors are influenced by the under- 
lying topological properties of the networks [EH]- Among 
many dynamical processes, a random walk on networks 
is fundamental to many branches of science and engi- 
neering, and has been the focus of considerable atten- 
tion g [H [ll, [H, [H, [13, El. As a fundamental dynam- 
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ical process, the random walk is related to various other 
dynamics such as transport in media [l6j , disease spread- 
ing [TtJ , target search [l8[ , and so on. On the other hand, 
the random walk is useful for the study of network struc- 
ture, in particular for the average path length [§, fill ]. 
It is thus of theoretical and practical interest to study 
a random walk on complex networks, revealing how the 
structure (e.g. structural heterogeneity) effects the dif- 
fusive behavior of the random walk. 

In addition to the scale-free behavior, it has also been 
observed that real networks display ubiquitous degree 
correlations among nodes [p| [2(J , which are usually 
measured by two quantities, i.e., average degree of near- 
est neighbors of nodes with a given degree [20] and Pear- 
son correlation coefficient [2l|, both of which are equiv- 
alent to each other. Degree correlations are important 
in characterizing network topology, according to which 
one can classify complex networks into categories [21|: 
assortative networks, disassortative networks, and uncor- 
related networks. For example, social networks are usu- 
ally assortative, while technological and biological net- 
works disassortative. Furthermore, degree correlations 
significantly influence the collective dynamical behav- 
iors, including intentional attacks on hub nodes [221. [23j]. 
games [24j . and synchronization [251 ] . to name but a few. 

In view of the importance of both the inhomogeneous 
degree distribution and degree correlations, some funda- 
mental questions rise naturally: In heterogenous corre- 
lated networks, how does the structural heterogeneity, 
characterized by the alterable degree exponent 7, affect 
the scaling character of the average path length d(N)l 
Does the relation between d(N) and 7 in uncorrelated 
networks also hold for networks with degree correlations? 
Is the behavior of random walks related to structural het- 
erogeneity in correlated networks? Such a series of im- 
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FIG. 1: (color online) The first four generations of the network 
construction for a special case of m = 2. 



portant questions still remain open. 

In this paper, we study the average path length of and 
a random walk on a family of deterministic treelike disas- 
sortative scale-free networks with changeable degree ex- 
ponent 7 G (2,1 + |^| ] . We choose deterministic net- 
works as investigation object, because they allow us to 
study analytically their topological properties and some 
dynamic processes running on them. Our exact results 
show that in contrast to the scaling obtained for uncorre- 
cted networks, in their full range of 7, the average path 
length of the considered deterministic networks grows 
logarithmically with the number of nodes, and that only 
partial scalings of the random walk depend on the degree 
exponent 7. 



II. THE RECURSIVE SCALE-FREE TREES 

In this section, we introduce a network model defined 
in a recursive way 12611 . which has attracted a great 
amount of attention [23, [H, El HH. We cal1 this 
model recursive scale-free trees (RSFTs) [26] ]. We inves- 
tigate RSFTs because of their intrinsic interest and be- 
cause these networks have general degree distribution ex- 
ponent 7 G (2,1 + . Moreover, RSFTs are determin- 
istic, which allows us to study analytically their topolog- 
ical properties and dynamical processes on them. They 
are therefore good test-beds and substrate networks. 

The recursive scale-free trees, denoted by R m ,t (jn is 
a positive integer) after t (t > 0) generation evolution, 
are constructed as follows. For t = 0, i? m ,o is an edge 
connecting two nodes. For t > 1, R m ,t is obtained from 
Rm,t-V For each of the existing edges in Rm,t—\i m new 
nodes are introduced and connected to either end of the 
edge. Figure Q] shows the construction process for the 
particular case of m = 2. 

According to the network construction, one can see 
that at each step ti (U > 1) the number of newly in- 
troduced nodes is L v (ti) = 2to(1 + 2m) ti ~ 1 . From this 



result, we can easily compute network order (i.e., the to- 
tal number of nodes) N t at step t: 

t 

N t = y £ j L v {t i ) = {2m+l) t + l, (1) 

ti=0 

Let ki(t) be the degree of a node i at time t, which 
entered the networks at step ti (U > 0). Then 

h{t) = {m+l) i - ti . (2) 

From Eq. (J2J), one can easily see that at each step the 
degree of a node increases m times, i.e., 

ki(t) = (m + l)ki(t-l). (3) 

RSFTs present some typical characteristics of real-life 
networks in nature and society, and their main topo- 
logical properties are controlled by the parameter m. 
They have a power law degree distribution with expo- 
nent 7 = 1+ 1 ' 1 ( 2m + 1 ) belonging to the interval between 
2 and 3 [H, H|. The diameter of RSFTs, defined as 
the longest shortest distance between any pair of nodes, 
increases logarithmically with network order [30], that 
is to say, RSFTs are small-world. The betweenness dis- 
tribution exhibits a power-law behavior with exponent 
7t = 2 (27[. In addition, RSFTs are disassortative, the 
average degree of nearest neighbors, k nD (k), for nodes 
with degree k is approximately a power-law function of 
k with a negative exponent [301 ] . 

After introducing the RSFTs, in what follows we will 
study the average path of the recursive scale-free trees 
and random walks on them. We will show that the expo- 
nent 7 of degree distribution has no qualitative effect on 
APL and mean first-passage time (FPT) for all nodes, 
but has essential influence on FPT for old nodes when 
the networks grow. 

III. AVERAGE PATH LENGTH 

We now study analytically the APL dt of the recursive 
scale-free trees R m ,t by using a method similar to but 
different from that proposed in Ref. 32]. It follows that 

dt = N t (N t - l)/2 ' (4) 

where $t is the total distance between all couples of 
nodes, i.e., 

in which rfy is the shortest distance between node i and 
j- 

Notice that in addition to the recursive construction, 
RSFTs can be alternatively created in another method. 
Given a generation t, R m ^+\ may be obtained by joining 
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In order to determine A*' 2 and A 2 ' m+2 , we define 



FIG. 2: (color online) Second construction method of RSFTs. 
The graphs after t + 1 construction steps, Rm,t+i, may be ob- 
tained by the juxtaposition of 2m + 1 copies of R m ,t, denoted 
by Rm t t (0 = 1, 2, 3, • • • , 2m, 2m + 1), which are connected to 
one another at the border nodes X (a) and Y (■). 



at hub nodes 2m + 1 copies of R m ,t, see Fig. [2] The sec- 
ond construction method highlights the self-similarity of 
RSFTs, which allows us to address dt analytically. From 
the obvious self-similar structure, it is easy to see that 
the total distance $t+i satisfies the recursion relation 



$t+i = (2m + l)$ t + A t 



(6) 



where A 4 is the sum over all shortest paths whose end- 
points are not in the same t branch. The solution of 
Eq. © is 



(2m + I)*" 1 *! + ^(2m + l)*~ r_1 A T • 



(7) 



The paths contributing to A t must all go through at least 
either of the two border nodes (i.e., X and Y in Fig. ^ 
where the different R B m t branches are connected. The 
analytical expression for A t , called the crossing paths, is 
found below. 

Let A" l/3 be the sum of all shortest paths with end- 
points in R^ t and R m t . According to whether or not 
two branches are adjacent, we split the crossing paths 
A"'^ into two classes: If R^t and R^.t meet at a border 
node, A"'' 3 rules out the paths where either endpoint is 
that shared border node. For example, each path con- 
tributing to A t ' 2 should not end at node X. If t and 

RP m t do not meet, A"'^ excludes the paths where either 
endpoint is X or Y. For instance, each path contribu- 
tive to A 2 ' m+2 should not end at nodes X or Y. We 
can easily compute that the numbers of the two types 
of crossing paths are m 2 + m and m 2 , respectively. On 
the other hand, any two crossing paths belonging to the 
same class have identical length. Thus, the total sum A t 
is given by 



At = (m 2 + m) A t 1,2 + ' 



; A 



2,m+2 



(8) 



0t = d lX 



(9) 



Considering the self-similar network structure, we can 
easily know that at time t+1, the quantity at+i evolves 
recursively as 

cr t+ i = (m+l)a t + m[a t + (N t -l)] 

= (2m + l)a t +m(2m+ 1)*. (10) 



Using do = 1 ; we have 

er t = (mt + 2m + 1) (2m + l)*" 1 . 



fill 



Having obtained at , the next step is to compute the quan- 
tities A]' 2 and A 2 ' m+2 given by 



At 1 ' 2 = d » 



[dix+djx) 
(N t — 1) £ d iX + (N t -l) d jx 

i^X j^X 

2(N t -l)a t , (12) 



and 



A 



2,m+2 



} J [dix + dxY + d jY ) 



ieR^i^X 

3^Rzi 2 ,3^y 

2(jV t - 1) a t + {N t - l) 2 , (13) 



where dxy = 1 has been used. Substituting Eqs. {T 
and (|13p into Eq. ©, we obtain 



At = 2m(2m+ l)(N t - l)a t + m 2 (N t - l) 2 

= m{2mt + 5m + 2)(2m + l) 2t . (14) 

Inserting Eqs. (fl"4]) for A T into Eq. ([7]), and using $1 = 
5m 2 + 4m + 1 , we have 

(2m + lY^ 1 

$ t = ^ '- [l + m + (2mt + 3m + 1) (2m + 1)*] . 

2 (15) 
Plugging Eq. (fl~5|) into Eq. ((4]), one can obtain the ana- 
lytical expression for dt- 

1 + m + (2mt + 3m + \)(2m + 1)* 
t_ (2m+l)[(2m+l) t + l] ' ^ 
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— O - Analytical Results for m = 3 






♦ Numerical Results for m = 3 
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FIG. 3: (color online) Average path length d t versus network 
order Nt on a semi-log scale. The solid lines are guides to the 
eyes. 

which approximates 2m+i * n ^ ne mnn ite t, implying that 
the APL shows a logarithmic scaling with network order. 
Therefore, RSFTs exhibit a small-world behavior. No- 
tice that this scaling has been seen previously in some 
other deterministic disassortative scale-free networks in 
the same exponent r ang e, such as the pseudofractal scale- 
free web studied in [33|, [34| and the "transfractal" recur- 
sive networks addressed in [35[ . 

We have checked our analytic result against numerical 
calculations for different m and various t. In all the cases 
we obtain a complete agreement between our theoretical 
formula and the results of numerical investigation, see 
Fig.H 

The logarithmic scaling of APL with network order in 
full rage of degree exponent 7 shows that previous rela- 
tion between APL and 7 obtained for uncorrelated scale- 
free networks 0,[1| is not valid for disassortative scale- free 
networks, at least for RSFTs and some other determinis- 
tic scale-free graphs. This leads us to the conclusion that 
degree exponent 7 itself does not suffice to characterize 
the average path length of scale- free networks. 

IV. RANDOM WALKS 

This section considers simple random walks on RSFTs 
defined by a walker such that at each step the walker, lo- 
cated on a given node, moves to any of its nearest neigh- 
bors with equal probability. 

A. Scaling efficiency 

We follow the concept of scaling efficiency introduced 
in [ll[ ■ Denote by the mean first-passage time (FPT) 
between two nodes i and j. Let Tu be the mean time for 
a walker returning to a node i for the first time after the 



walker have left it. When the network order grows from 
N to gN, one expects that in the infinite limit of N 

T^gN) ~ g^T^N), (17) 

where 8ij is defined as the scaling efficiency exponent. 
An analogous relation for Tu defines exponent 5u. 

One can confine scaling efficiency in the nodes already 
existing in the networks before growth. Let T^(gN) be 
the mean first-passage time in the networks under consid- 
eration, averaged over the original class of nodes (before 
growth). Then the restricted scaling efficiency exponent, 
Ay , is defined by relation 

VjigN) ~ g^T^N). (18) 

Similarly, we can define Xa. 

After introducing the concepts, in the following we will 
investigate random walks on RSFTs following a similar 
method used in [ll[. It should be mentioned that our 
motivation is different from that of [11]. In the work 
of [ll|, the authors analyzed scale- free networks with a 
single degree distribution exponent 7=1 + the pur- 
pose of that work is to find what is special about random 
walks on scale-free networks, compared to other types of 
graphs. Here we study random walks on scale- free net- 
works (RSFTs) with general 7 € (2, 1 + j|§] . Our aim is 
to study the effect of degree exponent 7 on random walks 
characterized by the scaling efficiency proposed in [ll[ . 

B. First- passage time for old nodes 

Consider an arbitrary node i in the RSFTs R m ,t after t 
generation evolution. From Eq. ([2]), we know that upon 
growth of RSFTs to generation t + 1, the degree fe« of 
node i increases to times, namely, from ki to (to + l)ki. 
Let the FPT for going from node i to any of the fc, old 
neighbors be T. Let the FPT for going from any of the 
mki new neighbors to one of the fcj old neighbors be A. 
Then we can establish the following equations 

{ T = T^fi + T^FlC 1 + A )> (10) 
\ A = 1 + T, 1 > 

which leads to T = 2m + 1. Therefore, the passage time 
from any node i (i € R m ,t) to any node j (j € R m ,t) 
increases 2m times, on average, upon growth of the net- 
works to generation t + 1, i.e., 

T! j (N t+1 ) = (2m + l)T ij (N t ). (20) 

Since the network order approximatively grows by 2to 
times in the large t limit, see Eq. (fTJ). This indicates that 
the scaling efficiency exponent for old nodes is Ay = 1, 
which is a constant independent of the degree exponent 

7- 

Next we continue to consider the return FPT to node i. 
Denote by the FPT for returning to node i in R m ^+i- 
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Denote by the FPT from j — an old neighbor of i 
(i G R m ,t) — to i, in R m j+i- Analogously, denote by 
Ta the FPT for returning to i in i? m ,t, and T yl the FPT 
from the same neighbor j, to i, in R m ,t- For R m ,t, we 
have 



T - — 



1 



E a 

jefii(*) 



1 



E r - 



(21) 



where f2j(i) is the set of neighbors of node i, which belong 
to R m .t- On the other hand, For R m j+i, 



T'~ = 



m + 1 



x 2 + 



+ 1 Kj ' J 



The first term on the right-hand side (rhs) of Eq. (|2"2"|) ac- 
counts for the process in which the walker moves from 
node i to the new neighbors and back. Since among all 
neighbors of node i, -^fr of them are new, which is obvi- 
ous from Eq. © , such a process occurs with a probability 
— ^tt and takes two time steps. The second term on the 

m+l " 

rhs interprets the process where the walker steps from i 
to one of the old neighbors previously existing in R m j 
and back, this process happens with the complimentary 
probability 1 - ^ = 

Using Eq. (|2"P)) to simplify Eq. ([2"2")l . we can obtain 

= ^±l Tll = ( 2m+ l)[i-M"»+i)/M2m+i)] Tii . (23) 
m + l 

In other words, 

T> l (N t+1 ) = (2m+l) x »T u (N t ), (24) 



, ln(ro+l) 
1 lu(2m+l) 



where the scaling efficiency exponent Ai 
1 — -^j- is an increasing function of 7. Thus, the more 
heterogeneous the network structure, the more easily for 
the walker to return to the origin when the networks grow 
in size. 



go to node j'; with the complimentary probability ^7^, 
the walker chooses uniformly a neighbor node except j', 
and spends on average time B in returning to i, then 
takes time T\ to arrives at node j'. 

In order to close Eqs. (|25|) and (j26|) . we express the 
FPT to return to i as 



1 „ h - 



T u (N t ) 

Eliminating T\ and R, we obtain 



1 



-B. 



T rr {N t )=kiT u {N t ). 
Combining Eqs. ©, ^ and we have 
r i / i /(JV t ) = 2(2m+l) t ~2JV t . 



(27) 



(28) 



Iterating Eqs. (|2"3"|) and 

are L^(e) (0 < e < i) nodes with T n - (m+1)t 
average of T,j over all nodes in R m ,t leads to 



we have that in R m ^ there 

2(2m+l) t 



(29) 

ffiere 
Taking 



1 



(2m + 1)* + 1 



4m(2m+l)* 4m(2m+l) f 



(m + l)* (m + l)'" 1 

(2m + l)*(m+ l) 4 - 1" 



4m + 4 
2m + 3 



(2m + 1)' 



m(2m + 3) 
4m + 4 



2m + 3 



(30) 



Equation (|30[) implies that 8a — 1, which is uncorrclatcd 
with the degree exponent 7. 

We continue to calculate in i? mi t, which is FPT 
from an arbitrary node i to another node j. Since each of 
the newly created nodes has a degree of 1 and is linked to 
an old node, the FPT Ti'j from node i' — a new neighbor 
of the old node i — to j equals Ty plus one and thus has 
little effect on the scaling when network order N is very 
large. Therefore, we need only to consider FPT X^-/ from 
i to j' — a new neighbor of j, which can be expressed as 



C. First-passage time for all nodes 

We now compute Tjiji , the FPT to return to a new 
node f £ R m ,t that is a neighbor of node i € R m ^-i- 
Denote by T Y the FPT from i to j', and B the FPT to 
return to i (starting off from i) without ever visiting j'. 
Then we have 



and 



Tyy =1 + 71, 



r, = I + *L-i(B + T 1 ). 



hi ki 



(25) 



(26) 



Equation (1261) can be interpreted as follows: With prob- 
ability £7 (fcj being the degree of node i in R m ,t), the 
walker starting from node i would take one time step to 



Notice that 



Substituting Eq. 
results in 



T -i 



T-i -i 



Tu + T. 



33' ■ 



l + T. 



33' ■ 



(31) 



(32) 



and Eq. 



for Tjijt into Eq. §1 



Ty + 2(2m+l) t -l^iV t , (33) 



where Eq. ([20]) has been used. Therefore, we have 



(Tij)t 



(34) 



which shows that mean transit time between arbitrary 
pair of nodes is proportional to network order. Equa- 
tion (|3"4")1 also reveals that 5y is a constant 1, which does 
not depend on 7. 
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V. CONCLUSIONS 

To explore the effect of structural heterogeneity 
on the scalings of average path length and random 
walks occurring on disassortative scale-free networks, we 
have studied analytically a class deterministic scale-free 
networks — recursive scale-free trees (RSFTs) — with vari- 
ous degree exponent 7. In addition to scale- free distribu- 
tion, RSFTs also reproduce some other remarkable prop- 
erties of many natural and man-made networks: small 
average path length, power-law distribution of between- 
ness distribution, and negative degree correlations. They 
can thus mimic some real systems to some extent. 

With the help of recursion relations derived from the 
self-similar structure, we have obtained the solution of 
average path length for RSFTs. In contrast to the well- 
known result that for uncorrelated scale-free networks 
with network order N and degree exponent 2 < 7 < 3, 
their average path length d(N) behaves as a double 
logarithmic scaling, d(N) ~ In In TV, our rigorous solu- 
tion shows that the APL of RSFTs behaves as a log- 
arithmic scaling, in despite of their degree exponent 
7 G (2,1 + ] . Therefore, degree correlations have a 
profound impact on the average path length of scale- free 
networks. 

We have also investigated analytically random walks 
on RSFTs. We have shown that for the full range of 
7, the mean transit time Tij(N) between two nodes av- 
eraged over all node pairs grows linearly with network 
order N. The same scaling holds for the FPT Tu(N) for 
returning back to the origin i after the walker has started 
from i. Thus, despite different 7, all the RSFTs exhibit 
identical scalings of FPT and return FPT for all nodes. 



On the other hand, for those nodes already existing in the 
networks before growth, the restricted scaling efficiency 
exponents are Xij = 1 and Xu = 1 — — zj, where Xij is 
not pertinent to 7, but Xu depends on 7. 

We should stress that our conclusions were drawn only 
from a particular type of deterministic treelike disassor- 
tative networks. It is still unknown whether the conclu- 
sions are also valid for stochastic disassortative networks, 
especially for networks in the presence of loops. But our 
results may provide some insights into random walk prob- 
lem on complex networks, in particular on trees. More 
recently, the so-called border tree motifs have been shown 
to be significantly present in real networks (36j . looking 
from this angle, our work may also shed light on some 
real-world systems. Finally, we believe that our analyt- 
ical techniques could be helpful for computing average 
path length of and transit time for random walks on other 
deterministic media. Moreover, since exact solutions can 
serve for a guide to and a test of approximate solutions 
or numerical simulations, we also believe that our vigor- 
ous closed-form solutions can prompt related studies of 
random networks. 
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